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The guiding lines of the lattice investigations on colour confinement are reviewed, to- 
gether with recent results. 

1. Introduction 

Colour is confined in nature. The experimental upper limit on abundance of quarks n q , 
on abundance of protons n p is 

^ < 10~ 27 (1) 

n p 

corresponding to Millikan like analysis of ~ lg of matter. If there were no confinement a 
conservative expectation in the cosmological standard model would be|TJ 

^ ~ 10- 12 (2) 
n p 

A factor 10~ 15 is convincing evidence for confinement. Its smallness also indicates that 
confinement cannot be explained by fine tuning of a small parameter, but only in terms 
of symmetry. An analogous situation exists in superconductivity 0. The persistence of 
electric currents during years in a superconductor is explained in terms of the change of 
symmetry of the ground state produced by the Higgs phenomenon. The lifetime of the 
current is determined by the characteristic correlation time of macroscopic fluctuations. 

Confinement must be explained in terms of symmetry and if QCD is the theory of 
hadrons it should have the relevant symmetry built in. 



2. Lattice QCD 

QCD most likely exists as a field theory due to asymptotic freedom, in contrast to other 
field theories like QED, which can only be viewed as effective models, and loose their 
validity at short distances. In rigorous terms this means that the euclidean Feynman 
integral which defines the theory 

27V 

exp(-pS) (3) 
admits a thermodynamical limit. 
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Z is a functional integral, and is defined as the limit of a sequence of discretized, 
ordinary integrals. Lattice is an approximant in this sequence, and a good approximation 
to the limit whenever the correlation length is large compared to lattice spacing, and 
small compared to the lattice size||. 

Physical quantities computed on the lattice, e.g. by numerical simulations, are therefore 
determined from first principles. Lattice formulation also defines the correct ground state 
(vacuum). The usual perturbative quantization in terms of Fock space of gluons and 
quarks is instead basically ill defined, and the instability of its ground state manifests itself 
by the existence of singularities (renormalons) in the resummation of the perturbative 
series @. Lattice is therefore specially suited to understand the structure of the theory. 

A deep issue in this understanding is the N c — > oo limit. The idea|5|] is that the 
limit N c — > oo of QCD, with A = g 2 N c fixed, defines a field theory, which contains the 
main physical features of QCD, including confinement of colour. An expansion in 1/N C , 
Nf I N 2 ... is well defined and is a correction to the limiting theory. A fermion loop in 
this expansion is roughly speaking 0(Nf /N^ 2 ) where V is the number of gluon-quark 
vertices. For V > 2 this is a small correction, and only the loops with V = 2 are 0(1). 
Quark loops are negligible except for the rescaling of A Lattt due to contribution of quarks 
loop with V = 2 to the (3 function. Numerical simulations on lattice do support this 
expectation. 

An alternative evidence in this direction is the solution of the U(l) problem|6|], where 
the mass of the rj', m v ' is related by l/N c arguments to the topological susceptibility x of 
the quenched QCD vacuum]?] 

2N fX = flm\ (4) 
where 

X = [ d i x(T(Q(x)Q(0))\0) (5) 



Q(x) is the topological charge density, Nf the number of light fermions, f n the pion decay 
constant. 

Eq.fJD predicts x(180MeV) 4 and this prediction is confirmed by lattice simulations ||. 
All this points to the independence of the confinement mechanism on N c , Nf. 
The change of symmetry leading to confinement should not depend on N c , Nf. 

3. Finite temperature QCD 

The standard way to define a static thermodynamics of QCD is to limit the euclidean 
time direction from to 

t E = l/T (6) 

with periodic boundary conditions for bosons, antiperiodic for fermions. 

In lattice QCD an asymmetric lattice is used Ng x N t , with iV s ^> N t . At a given value 
of the lattice unrenormalized coupling (3 = 2N c /g 2 the temperature will be defined by 
eq.(|B|) as 

<P)Nt = f (7) 
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in terms of the lattice spacing a{f3). 

The value of a{(3) in physical units, depends on (3 via renormalization group. At one 
loop 

a((3) = -^—exp(-b f3) (8) 

^Latt 

where bo > is the negative of the first coefficient of the perturbative expansion of the (3 
function. 

At a given (3 then 

T=^exp(6 /5) (9) 

Temperature increases exponentially with f3: the deconfining transition from hadrons to 
quark gluon plasma will take place at some value of j3. For T > T c , or in the weak coupling 
regime, the order parameter of the theory, the Polyakov line (P) is ^ 0, and quarks and 
gluons exist as particles ||. For T < T c , (P) = 0, and there is confinement. In general 
(P) is related to the chemical potential of a single quark, 7, by the equation 

(P) = exp(- 7 ) (10) 

so that the confined phase corresponds to 7 = 00: an infinite amount of energy is needed 
to create an isolated quark. From the point of view of symmetry T < T c is the disordered 
phase, or the strong coupling regime. 

If confinement is due to symmetry, the obvious question is: what is the symmetry of a 
disordered phase? As we shall see the key word in this respect is "duality" . 

Before going to duality we underline two points: 

a) (P) is an order parameter only in the absence of quarks. In view of the N c — > 00 
argument presented above, such situation is not satisfactory. A order parameter 
should exist for the two situations. 

b) Evidence from lattice exists that confinement takes place in QCD. Below T c the QQ 
potential V(R), as determined from Wilson loop, is [I 



V(R) = aR a = string tension (11) 

Above T c the force is a screened Coulomb force. Moreover in the confined phase the 
linear potential is related to the existence of chromoelectric flux tubes joining the 



QQ pair, whose energy is proportional to the length |Tl|J12 
4. Duality 

Duality is a typical property of (d + 1) dimensional systems which can have spatial (d 
dimensional) configurations with non-trivial topology and conserved topological charge. 
At low g 2 the system is described in terms of local fields $, the symmetry is identified in 
terms of order parameters ($) and there is order. At high g 2 , ($) — > 0, and disorder sets 
in. 
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However a dual description exists, in which the extended structures with topology are 
described by local fields, /i, the original $ fields become non local, and the effective 
coupling constant is g' ~ l/g 2 - In this description the original disordered phase is ordered 
and viceversa. The dual symmetry is described by (/j,) which is called a disorder parameter. 

The prototype model of duality is the 2 dimensional Ising model, where the topological 
configurations are kinks | 



Other examples are liquid He^, with vorticesJT^], the Heisenberg ferromagnet ||16|| , with 
non abelian vortices or Weiss domains, and the usual 3+1 dimensional gauge theories, 
with monopoles. U{1) is an example [jlTf, SUSY QCD[T^] is another and finally QCD, as 



we shall show in what follows. 

The strategy that we adopt is to identify the dual symmetry, and to construct the 
disorder parameter (/i) in terms of the original local fields $, which in QCD are gluons 
and quarks. This strategy is inspired by ref. |Ll| and has been adapted on lattice by our 



group. An alternative strategy consist in performing the transformation to dual [0,0] , 
but is less convenient in numerical investigations. 

5. Monopoles in QCD 

A guess of the dual symmetry in QCD is that the vacuum in the confined phase is a dual 
superconductor |20 . The phenomenological motivation is that in a dual superconductor 



dual Meissner effect will constrain the chromoelectric field acting between a QQ pair into 
an Abrikosov like flux tube, with energy proportional to the length, or to the distance 

V{R) = oR (12) 

Monopoles are objects with nontrivial topology, when described in terms of the gauge 
fields. The signal of dual superconductivity should be the nonvanishing of the vacuum 
expectation value of an operator carrying magnetic charge, the analog of the Ginsburg 
Landau charged Higgs field in ordinary sup er conduct ivitypT] . 



Monopoles in QCD are well understood [^,0]. They can be exposed by a procedure 
known as abelian projection|24]]. We shall sketch it for SU(2) gauge fields: the extension 
to SU(N) is trivial. 

To each operator <3>(x) = $(x) • | in the adjoint representation a gauge transformation 
R(x) can be associated, which brings $ parallel to 3 axis, or $(x) to a diagonal 2x2 
matrix, proportional to cr 3 . R(x) is called "abelian projection" and is defined everywhere 
except at zeros of $(x). We shall define for convenience 

$(x) = iH- (13) 

\$(x) 



The field strength || 



F IIV = Q- - -$ ■ (D^ A £>„$) (14) 

with G^y the usual field strength and = — gA^A the covariant derivative, is a 
gauge invariant, colour singlet operator. By construction the bilinear term in A^ cancels 



between the two terms in eq.(|i~4[) and in the abelian projected gauge 

Fpv = d„Al - d v A\ (15) 
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looks like an abelian field. Its dual F* v 



F pu defines a magnetic current jff as 



J n 



d v F, 



which is identically conserved 
*rft = 



(16) 



(17) 



A magnetic, gauge invariant U(l) symmetry is thus associated to each operator $(a;) 
belonging to the adjoint representation. 

jjf is not trivially zero: monopoles charges can exist at the sites where <&(x) = 0, where 



the abelian projection R(x) is singular [23 1 . 

Dual superconductivity will occur when magnetic charges belonging to the above U(l) 
condense in the vacuum, and the corresponding electric field, i.e. the chromoelectric field 
parallel to the 3 color axis after abelian projection, will experience dual Meissner effect. 

The disorder parameter (p) describing monopole condensation can be constructed with 
a technique tested in a variety of known system fT5HT7j|, following the strategy outlined 
above. We will not go into the details of the construction, for which we refer to|25|. 

The transition to superconductor is signalled by a sharp negative peak in the plot of 
the quantity p = -^ln(p) versus (3 (or vs T). The peak becomes sharper as the space 
volume of the system becomes larger. A typical behaviour is shown in fig.(l). 
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Figure 1: (p) vs. (3 for SU(2) gauge theory. 
Plaquette projection, lattice 16 3 x 4. 

A finite size scaling analysis allows to extract the critical indices and the decoupling 
, lq- 17|| . Indeed, for dimensional reasons 



temperature [ 



L d <$>(- 

A Ng' 



;i8) 



where a is the lattice spacing, A is the correlation length, N s the extension of the system. 
At p c , X diverges as A ~ (j3 c — f3)~ u , with v the correlation critical index, and a/A C 1. 
In the limit a/ A — >• the scaling law follows for p = ^ln(/i) 



P 

-Xjv 



$(L^(/3 c -/3)) 



(19) 



6 




Figure 2. Quality of the scaling for SU(2) Figure 3. Quality of the scaling for SU(3) 
gauge theory. Plaquette projection, N t = 4. gauge theory. Polyakov projection, N t = 4. 



The quality of scaling is shown in fig.2 for SU{2) and in fig. 3 for SU(3). Scaling exists 
only for the proper values of /3 C , u, which can thus be determined. 

These quantities can also be determined by different methods0. They agree within 
errors. For SU (2) we find 

v = 0.63(5) 
P c = 2.30(2) 
5 = 0.20(8) 

For SU(3) 

v = 0.33(2) 
P c = 5.70(3) 
5 = 0.54(4) 

6. Discussion: open problems. 

There is definite evidence that confinement is produced by dual superconductivity. 
What has not been specified above is the field $(x). In principle monopoles exist for any 
choice of it: a functional infinity of monopoles species. There is no argument a priori that 
all of them should condense. In ref|24]] t'Hooft guesses that all monopoles are physically 
equivalent, and all of them should condense in the confined phase. What we observe on 
lattice supports this view. A few choices for $ that we have tried look perfectly equivalent. 

An alternative attitude exists in the literature, and is that the so called "maximal 
abelian projection" defines monopoles which are more equal than others. In fact with 
that choice the abelian field dominates, since after the projection all the fields are oriented 
in its direction to 90%. This choice can be more convenient than others if one looks for 



N t = A 



N t = 4 , 
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effective lagrangians. In principle this convenience does not put any limitation to the 
symmetry pattern of confinement. 

Most probably a better and more syntetic understanding of this symmetry is needed. 
There must be a reason why all abelian projections are equivalent. However there is no 
doubt that dual superconductivity is the mechanism of colour confinement. 

The part of the work reported due to our group has been done in collaboration with 
L. Del Debbio, G. Paffuti, P. Pieri, B. Lucini, D. Martelli in the last few years. Their 
contribution was determinant to the results. 
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